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Abstract 

We compare two different approaches to investigations of many-electron systems. 
The first is the Hartree-Fock (HF) method and the second is the Density Functional 
Theory (DFT). Overview of the main features and peculiar properties of the HF 
method are presented. A way to realize the HF method within the Kohn-Sham (KS) 
approach of the DFT is discussed. We show that this is impossible without including 
a specific correlation energy, which is defined by the difference between the sum of 
the kinetic and exchange energies of a system considered within KS and HF, respec- 
tively. It is the nonlocal exchange potential entering the HF equations that generates 
this correlation energy. We show that the total correlation energy of a finite elec- 
tron system, which has to include this correlation energy, cannot be obtained from 
considerations of uniform electron systems. The single-particle excitation spectrum 
of many-electron systems is related to the eigenvalues of the corresponding KS equa- 
tions. We demonstrate that this spectrum does not coincide in general with the 
eigenvalues of KS or HF equations. 



PACS number(s): 31.15Ew, 31.50.+W 



Typeset using REVT[t;X 



*E-mail: vrshag@thd.pnpi.spb.ru 



1 



I. INTRODUCTION 



The Hartree-Fock method (HF) originate from the pioneer works of Hartree [1] and Fock [2]. 
Since then an impressive progress has been achieved by using HF as a basis when calculating the 
ground state energy and constructing the Random Phase Approximation (RPA) and the Random 
Phase Approximation with Exchange (RPAE) to investigate the dynamic properties of multi-electron 
objects, such as atoms [3], molecules, clusters and fuUerenes. Unfortunately, HF method has a number 
of difficulties in its application. For instance, one needs a whole set of single-particle wave functions 
to calculate the single-electron nonlocal potential. It is also necessary to follow a rather intricate 
procedure to include the correlation corrections that are beyond the HF framework. As a result, 
in the case of complex many-electron systems, starting with molecules, that include at least several 
atoms, the calculations become too comphcated. 

The Density Functional Theory (DFT) that was invented in the early 1960s has, generally speak- 
ing, overcome these difficulties. DFT has reduced the problem of calculating the ground state 
characteristics of a many-electron system in a local external field to the solution of the Hartree-like 
one-electron equations [4-6]. At first DFT was limited to consideration of only the ground states 
of many-electron systems, leaving aside their dynamic properties, which are closely related to the 
systems behavior in the time-dependent external fields. This limitation of DFT was surmounted by 
Runge and Gross [7] by transforming DFT into the so-called Time- Dependent Density Functional 
Theory (TDDFT), which in fact is a quite natural generalization of DFT. 

Both DFT and TDDFT are based on the one-to-one correspondence between the particle densities 
of the considered systems and external potentials acting upon these particles. For definiteness, we 
assume that the time dependent part of the considered electronic system's density p{r,t), which 
is created under the action of an external time-dependent field XvextiT^,t), is developing from the 
system's ground state. Consequently, one can conclude that the wave function and each observable, 
describing the system, can be presented as a functional of the density. Examples are the functionals 
of the ground state energy E[p] and of the action A[p] [5-7]. 

DFT and TDDFT estabhsh an exact correspondence between a system of many interacting par- 
ticles and a fictitious non-interacting Kohn-Sham (KS) system [5-8]. The main advantage of this 
"mapping" is that KS system is described by a set of Hartree- type single particle equations, which 
is not too difficult to solve. As a result, after solving these equations in the case of an external 
time-independent potential, one can obtain and predict, at least in principle, the atomic, molecular, 
cluster and solid bodies binding energies, phonon spectra, activation barriers etc. (see e.g. [6,8]). The 
same is true in the case of a time-dependent potential when the solution yields the single-particle 
and collective excitation spectra (see e.g. [8-11]). 

Unfortunately, the one-to-one correspondence proved in [4,7] establishes the existence of the 
functionals only in principle, leaving aside a very important question on how one can construct them 
in reahty. This is why the success of DFT and TDDFT strongly depends upon the availability of 
good approximations for the functionals. A systematic way to construct the required functionals and 
the corresponding response functions and effective interaction was suggested in a number of papers 
[9-11]. 

Since the HF and the Kohn-Sham theory play an important role in the many-body physics, 
it deserves comprehensive examinations in order to reveal their shortcomings and advantages as 
compared to one another when applied to studies of the many-electron systems. Recently, a key 
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element of DFT, the action functional A[p\, was suspected to be ill-defined because of a number 
of contradictions in it. It appeared that the contradiction, which came from the analysis of the 
causahty and symmetry properties of the hnear response functions in TDDFT, plays the central role 
in creating other contradictions [8,12,13]. Fortunately, these contradictions were successfully resolved 

recently [11,14-16]. 

Exhaustive calculations of the ground state energies of different atoms within the HF approach 
and that based on the density functional corresponding to the HF approximation (DHF) have shown 
that there exists an inequality, 

Edhf > Eqep > Ehf, (1) 

where Edhf is the ground state energy calculated in DHF, using a KS model system with the exact 
same density as that of some given HF ground states [17-19]. Here Eqep denotes the energy from the 
optimized potential method (0PM) [17,20], and Ehf denotes the HF ground state energy (see e.g. 
[17-19]). One can assume that inequality Eqep > Ehf results from the fact that Eqep is calculated 
with the restrictions associated with local potentials while calculations of Ehf are in this sense 
unrestricted (see e.g. [21]). But such considerations could justify only the relation Eqep > Ehf- For 
instance, if one applies this consideration to the functional of the KS kinetic energy Tk[p] one would 
expect that Tk[p] > T^^[p]. Quite surprisingly, one has instead T^^^lp] > Tk[p] ( sec e.g. [17]). Here 
T^^lp] is the HF functional of the kinetic energy. Then, one could expect, however, that instead the 
following relation holds (see, e.g. [18]) 

Edhf — Eqep — Ehf, (2) 

since a universal Hohenberg-Kohn (HK) functional is defined for HF ground states because of the 

constrained-search derivation [22]. Note, that this procedure also does not define explicitly the 
HF functional F^^[p] whose minimization leads to the HF ground state. The obvious contradiction 
between inequality (1) and equality (2) casts doubts on the reliability of the KS scheme and probably 
points to the possible contradictions in DFT based on the suggestion that an exact local exchange 
potential does not exist for ground states of typical atoms (see [18,19] and references therein). To solve 
the problem, we have to clarify the relationships between the HF exchange and kinetic functionals 
and the corresponding exchange and kinetic functionals of the DFT. 

Without resorting to a numerical analysis, by clarifying the relationship between the nonlocal 
exchange HF potential and the local exchange KS potential, we show that inequahty (1) and equality 
(2) have to be replaced by the inequality 

Eqep > Edhf = Ehf- (3) 

We also consider the important problem of how to relate the single particle spectrum of a many 
electron system to the corresponding KS eigenvalues. This consideration enables us to clarify the 
relationship between the single particle spectra obtained in HF calculations and in the KS calcula- 
tions. 

The structure of the paper is the following. In Section II, for the reader's convenience we review 
the main features of the HF method. We list only the main peculiarities of the HF method, which 
are directly determined by the nonlocal HF potential. One cannot expect to reproduce these within 
the Kohn-Sham theory of DFT. The explanation of Eqs. (1), (2), and (3) presented in Section HI 
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consists in demonstrating that the contradictions come from omitting the correlation energy, which is 
defined by the difference between the sum of the kinetic and exchange energies of a system considered 
within DFT and HF approaches, respectively. Consequently, this difference, which has to be treated 
as a correlation energy F^^[p], gives rise to a specific correlation energy, because the nonlocal HF 
potential ensures the deeper minimum value of the ground state energy than the local potential 
of 0PM. We also clarify the relationships between HF exchange and kinetic functional and the 
corresponding exchange and kinetic functionals of the DFT. Our consideration also shows that this 
correlation energy F^^[p\, being the specific feature of a finite system, cannot be represented by the 
local density approximation based on calculations of a homogeneous electron system. Therefore, the 
correlation energy of a finite electron system, which has to include F^^[p], cannot be obtained from 
considerations of uniform electron systems, although this way to construct Ec[p] is generally used 
[6] . In Section IV, we consider the single-particle excitation spectra of many-electron systems and 
relate them to the eigenvalues of the corresponding KS equations. We show that the single-particle 
spectra do not coincide with the eigenvalues of the KS or HF equations. A possible inconsistency in 
the Kohn-Sham theory based on calculations of the single-particle spectra [18,19] is also clarified in 
this Section. Finally, Section V contains our concluding remarks. 

II. THE HARTREE-FOCK METHOD 

Hartree's equations, as Fock's generalization of them, were derived for any atom with two and 
more electrons. In fact, they are applicable to any many-electron system in a static external potential 
v{r). Assuming that the total electron wave function \E'(ri, ...rjv) can be approximated as a product 
of one-particle functions 0f(r), ^^(ri, ...r^v) — Il^=i(l>f{^i): Hartree suggested [1] the following 
equation for 0|^(r)^: 

0f (r) = 0f (r). (4) 

Here is the total number of electrons. For an atom one has v{r) = —Z/r, where Z is the nuclear 
charge. An electron in an i-state is under the action of the so-called self-consistent field, which is 
determined by all electrons but the ith one. Excluding the term j — i from the sum in (4) eliminates 
the action of the i-clcctron upon itself. In Hartree approximation the system's electron density p^(r) 
is determined by the relation 

p^(r) = ^n,|0f(r)p, 
j=i 

where rij are the occupation numbers. 

Hartree's approach has two obvious shortcomings: the total wave-function violates the exclusion 
principle and, since the self-consistent potential acting upon the ith electron depends upon i, the 



^Atomic system of units is used in this paper: e = m = h = 1, where e and m are electron charge and 
mass, respectively. 
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functions 0f^(r) are non-orthogonal. These shortcomings were ehminated in [2] by constructing the 
total wave function as an anti- symmetrized {A) product, 

N 

^^^(ri,...r^) = An0r(r.), 

which automatically satisfies the exclusion principle. This wave-function's choice leads to the follow- 
ing ground state energy Ehf 



E 



HF 



^ .7 = li=l 



^-(ri)0:(r2)0,(ri)0,(r2) 
|ri - rsl 



dridr2 



(5) 



0*(rO0:(r2)0.(rO0,(r2) 



ri - r2 



dvxdvi + Y^nj / 0*(r)0j(r)v(r)dr, 

7=1 



and HF system of equations: 



0. 



HF I 



(6) 



TV 

+ E 



r — r' 



HF iHF I 



These equations differ from an ordinary Schrodinger equation in two essential aspects: they are non- 
linear in 0f^^(r) and the second term under the sum that represents the so-called Fock's potential is 
non-local. 

The non-linearity of Hartree and HF equations leads to non-uniqueness of their solutions (see [23] 
and references therein): in some cases, a given energy value can be represented by more than one set 
of solutions which deliver this value with high numerical accuracy. Furthermore, the nonlocal HF 
potential leads to a number of very specific features of the ground state wave function, which cannot 
be reproduced by the KS equations, see Sec. HI. 

1. The asymptotic behavior in r of (f)^^{r), contrary to the case of an ordinary one-particle 
Schrodinger and also Hartree equations, is not determined by Ef^ and does not behave as 



exp(-J2|Ef^|r) 



(7) 



On the contrary, it was shown in [24] that at r ^ oo the function (f)f^{r) is not determined by E^^ 
but is of the form 



0r(r)|.^oo ~ E,^^[exp(-^/2^r)]/r'+^ 



(8) 



where Ep^ is the energy of the so-called Fermi-lcvel (with wave function (j)p^{r)), which is the 
smallest binding energy among all E^ in the considered system, and 



Ci-T. j ci>r*i^yYi^ir/r)<j>f^ir)dr. 



(9) 
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Here Y;m(r/r) is the order spherical polynomial. The alteration of the asymptotic behavior 
dramatically affects, for example, the probability of penetration of atomic electrons via a potential 
barrier created when a strong electric static field is applied to this atom [25]. 

2. The number of zeroes in the radial part of the function (f)^^(r) is not equal, according to 
[26] and contrary to the case of an ordinary one-particle Schrodinger and also Hartree equations, to 
the so-called radial quantum number n,. = n — / — 1, where n is the principal and / is the angular 
momentum quantum numbers. The extra zeroes are located at big distances and their presence can 
be explained qualitatively in the following way. At large distances (f)f^{r) is determined mainly by 
the admixture of (f)p^{r) to 4>f^{r). The latter function adds its remote zeroes to that of 0f (r). 

3. In the presence of the HF potential the electron velocity operator v, which in quantum 
mechanics is determined by the relation v = i[H^^ , r], is not equal to p (note, that m = 1), contrary 
to the case of a local potential. This drastically affects, for example, the photoionization cross section 
calculations: the amplitude of this process becomes different in the so-called length and velocity forms 
and the golden sum rule is violated [3] . It is of interest to note that this violation is very serious and 
can be eliminated only when the corresponding calculations are performed within RPAE framework. 

4. One-electron Green's function G^^(ri,r2) in HF has only one expression via (j)^^{r), namely 

Ge E,-E + i5 ' 

while in case of the spherical symmetry there exist two equivalent expressions for local potentials 

where r>(<) is the larger (smaller) value of ri,2 and XE(r>) is the solution of the same Schrodinger 
(and also Hartree) equation as </>£'(r<), but irregular at r = 0. Note, that (J)e{i^) is regular at r — > 0. 

5. Another specific feature of the HF potential is the behavior as a function of the incoming 
electron energy E of elastic scattering phase shifts Si{E) of a partial wave with angular momentum 
I. As is well known, if the wave function's phase of electron scattering upon a local potential v{r) 
is normalized in such a way that Si{E)\e^oo 0, the relation Si{0) — qitt holds, where qi is the 
number of the electron's vacant bound states with angular momentum / in this potential (see e.g. 
[27]). The situation for the HF potential is qualitatively different. With the same normalization 
Si^^{E)\e^oo another relation holds: 5/^^(0) = {qi + ti)7i, where ti is the number of occupied 
bound states with angular momentum I in the HF potential [3,23]. 

III. THE KS LOCAL EXCHANGE POTENTIAL AND THE HARTREE-FOCK 

NONLOCAL ONE 

Consider a many-body system composed of interacting electrons and moving in an external 
potential v{r), with Hamiltonians H and Hy that are of the form 

H^f + U, (12) 

and 
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(10) 



H^^f + U + V. (13) 
Here T, U, and V are given in terms of the field operators ■0(r), and T is the kinetic energy operator 

T--lJr{r)V'i^{r)dr, (14) 

U is the interelectron interaction energy operator 

i/ = i/^:(£i)»(^H)Mrf,,rf,, (15) 

2 J |ri-r2| 
and V is the potential energy operator in the field v{r) 

J i/;*{r)iP{r)v{r)dr. (16) 

Hohenberg and Kohn proved that there exists a one-to-one correspondence between v{r) and the 
single-electron density p(r), that is v{r) is a unique functional of the density, apart from an additive 
constant [4]. Because v{r) fixes the Hamiltonian Eq. (13), the ground state wave function is a 
functional of the density, and one can define the energy functional [4] 

E[p] = F[p] + I v{r)p{r)dr, (17) 

with F[p] being a universal functional independent of the number of particles and the external 
potential v{t). Note that the functional F[p] can also be obtained on the basis of the constrained- 
search derivation of Levy, when it is defined by the minimum over the set of all normalized A^-electron 
wave functions that determine the density p(r) [22] 

F[p] = min (^^jf + t/j^^) . (18) 

Provided the density is exact, p(r) = po(r), and normalized by the number of particles N, N — 
J p{r)dr, the functional E[p\ assumes its minimum value that is equal to the ground state energy Eq 

Eo = F[po] + J v{r)po{r)dr. (19) 

Here the ground state density Po{r) is determined by the equation 

SF[p] 



dp{r) 



+ v{r)^p, (20) 



and p is the so-called chemical potential. The number p ensures the conservation of the number of 
particles A^. It is clear that if U = 0, [see Eq. (13)], the functional F[p] reduces to the functional 
Tfc[p] of the kinetic energy of noninteracting particles. 

Equation (20) can be replaced by a set of Hartree-like single-particle equations [see Eq. (4)], 
which permits the calculation of the density p(r) [5] 
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A 



(21) 



with the potential Vl — Vh + Vxc + The density p(r) is given by 

p(r) =^ni|0,(r)|2. 

i 

Here Vff is the Hartree potential [see Eq. (4)], 

bp\v) J |r — ri| 

where is the Hartree functional 

1 /•p(ri)p(r2) 



(22) 



(23) 



2 7 I" — i"i 



(24) 



and the exchange-correlation potential, Vxcir) = 6Fxc/Sp{r), is obtained from the exchange- 
correlation functional, 



FM = F{p] - T,{p] - FhIp]. 



(25) 



In the same way, it is possible to define separately F^ and the correlation functional Fc using the 
equation, which can be considered as a definition of Fc [5] 



Fc[p] = F,c[p] - F,[p]. 
Here F^lp] is the exchange functional [10,28] 

1 f Xo(ri, r2, iw) + 2'Kp{ri)6{w)6{ri - ra) dw 



F.[p] 



I 



k.i 



|ri - V2\ 



|ri - r2| 

and xo is the linear response function, which is of the form 
Xo(ri,r2,a;) = ^ nj(l - nfc)(/)*(ri)0i(r2)0fc(r2)0fc(ri) 

i.k 



27r 
dridr2, 

1 



dridr2 



(26) 



(27) 



u -ujik + ir] 00 + ujik - ir) 



(28) 



with cuik defined as cuik — Sk — Si, and functions 0fc(ri) being respectively the one-particle energies 
and wave functions of Eq. (21); 77 — > 0. Having the functional F^ given by Eq. (27), we can construct 
the HF-hke scheme based on Eq. (21) with the local KS exchange potential [10,28] 



K(r) 



5 



5p{v) 



(29) 



The functional of density F^^ corresponding to the HF approximation can be defined as follows 
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Eq. (12) and yielding the density p(r). This determinant is composed of the single-particle wave 
functions, which are solutions of a single-particle equation like Eq. (21) provided that VL(r) is a local 
potential. We denote this determinant by ^^^[p\. On the other hand, the same wave function in the 
frame of the optimized effective potential [20] is determined by the local exchange potential Vqpm- 
Because of the one-to-one correspondence between this wave function and the local single-particle 
potential [4] we can conclude that Vqpm coincides with given by Eq. (29). This result is in accord 
with the direct comparison of with Vqpm that confirms their equality, V^(r) = VopM(r) [10,28]. 
Therefore, the ground state energy Eqep of a many-electron system calculated in 0PM has to be 
equal to the corresponding energy calculated with the functional given by Eq. (30), E^sx, that is 
Eksx — Eqep- 

Now we turn to the HF calculations based on the constrained-search formulation of DFT, which 
leads to the HF functional F^^\p\ [22]. As it was mentioned in Section II and in accord with 
the constrained-search formulation [22], the HF functional is determined by the minimum of the 
expectation value F^^[p] — {^[p]\H\^[p]) when the minimization runs over all single-electron Slater 
determinants yielding the density p(r). The determinant, which yields the minimum value, is denoted 
as '^hf[p\- This construction of the functional F^^[p] makes it clear that the ground state energy 
Edhf of an electron system calculated with F^^[p\ is equal to E^f- Note that the determinant 
'^hf[p\ is composed of the eigenfunctions of Eq. (6) incorporating the Fock nonlocal potential. As 
a result, we can at first define the HF functional of the kinetic energy 



We will now show that T^^[p] ^ Tfc[p], and our proof proceeds by reductio ad ahsurdum. Note 
that there is a one-to-one correspondence between p(r) , the single-particle potential and the ground 
state wave function [4]. Because Tk[p] is a unique functional of the density, each value of Tk[p] 
fixes the only ground state function. This function is a single Slater determinant yielding the 
density p(r) normalized to electrons and characterized by the single-particle wave functions 4>i and 
eigenvalues Sj, given by Eq. (21). Assume that T^^[p] — Tj.[p]. Then the HF wave functions ^^^[p] 
and the eigenvalues E^^ given by Eq. (6) must be equal to $'^'^[p] and the eigenvalues Si given by Eq. 
(21). As it was mentioned in Sec. II, Eq. (8), all of the HF single-particle functions have the same 
long-range behavior which is dominated by the smallest orbital energy [24], while, generally speaking, 
the eigenvalues E^^ are not degenerate. On the other hand, eigenfunctions of Eq. (21) can have 
such a behavior only if the eigenvalues £j are degenerate. In order to eliminate the contradiction we 
have to recognize that T^^[p] ^ 7fe[p], and some of densities, which are the ground state densities of 
single-particle Schrodinger equations with nonlocal potentials are not noninteracting v representable 
in the sense that they cannot be the ground state density of Eq. (21) with some local potential f (r). 
Note that examples of densities which are not v representable were given in Ref. [29]. Obviously, a 




(31) 



and then the HF functional of the exchange energy 




(32) 
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one-to-one correspondence between nonlocal potentials and local ones does not exist [30]. Therefore, 
if we would even have the functional F^^[p] explicitly we would not have been able to calculate the 
HF ground state within the KS theory because the HF ground state cannot be reproduced by Eq. 
(21) with a local single-particle potential. On the other hand, the HF ground state can be obtained 
by solving the Hohenberg-Kohn equation (20) with F[p] = F^^[p]. These facts do not imply any 
contradictions. The Kohn-Sham theory and Eq. (20) lead to the same result if the density p(r) 
of system under consideration is v represent able. Otherwise one has to use Eq. (20). A detailed 
investigation of this important problem will be published elsewhere. 

The functional T^^ given by Eq. (31) is determined in the domain of the HF densities phf 
generated by the non-local Fock potential entering Eq. (6), while is determined in the domain of 
the densities pi given by Eq. (21). Note that both of the densities are normalized to same number of 
electrons. Approximating p^F by p^, one has to obtain < T^^ because of the minimum principle 
which ensures that the Slater determinant delivers the minimum value of the kinetic energy 
Tjfc. This result is confirmed by numerical calculations [17]. On the other hand, due to the minimum 
principle F^^ < F^^. Following Ref. [17], we can define the HF exchange-correlation functional 

Fx7 - -n-F^^ Tr + F^ - r„ (33) 
and introduce the HF correlation functional F^^ 

Fr - F[^^ - = Tr + Fr -n- f^. m 

We note that Eqs. (33) and (34) cannot be considered as exact definitions of F^^ and F^^ because it 
is impossible to approximate exactly Phf^Y Pl- It follows from Eq. (34) that the correlation energy 
F^ exists. This energy must be taken into account when constructing constructing approximations 
of F^^ within the KS framework. Omitting F^ ■, one obtains Eq. (1), while properly constructing 
F^^[p\ based on the constrained-search formulation of DFT [22], one obtains Eq. (3) rather then 
Eq. (2). 

In the case of He atoms the HF potential acting on the occupied states is a local potential. 
Therefore, the kinetic energy functional T^Ip] coincides with Tk[p\ leading to the equahty F^Ip] — 
Fx^[f^-i Eqep = EuHF = Ehf in line with Eqs. (2) and (3). This observation is confirmed by 
numerical calculations [18]. For all other atoms, the HF potential is nonlocal, and F^^[p] contributes 
to the ground state energy. Obviously, one has Eqep > Edhf = Ehf in accordance with Eq. (3) 
and the numerical calculations [18]. In the case of a homogeneous system, both single-particle wave 
functions, 0f^^(r) given by Eq. (6) and the ones 0i(r) given by Eq. (21), are plane waves. As a result, 
the energy Eksx — Ehf (though the Fock potential, Eq. (6), is a nonlocal) and the correlation energy 
Ec = E — Eksx = E — Ehf is uniquely defined. This consideration shows that F^^[p], being the 
specific feature of a finite system, cannot be approximated by the local density approximation based 
on calculations of a homogeneous electron system. Thus, Ec[p] of a finite electron system, which has 
to include F^lp], cannot be obtained from considerations of uniform electron systems, although 
this way to construct Ec[p] is generally used [6]. As a result, we are led to the conclusion that the HF 
method has unique peculiarities which can be considered as useful rather then as shortcomings when 
investigating many-electron systems. Moreover, it is of great importance to check in experiments 
whether the wave function of a real physical system has the behavior resembling the behavior of 
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the HF wave function given by Eq. (8) because this point is closely related to the problem of the 
f-representability. As mentioned in Sec. II, this HF-like behavior can be checked measuring the 
probability of penetration of atomic electrons via a potential barrier [25] . 



IV. SINGLE PARTICLE SPECTRA 

In order to calculate the single-particle energies e,j of a many-electron system we have to con- 
sider the ground state energy as a functional of the density p and of the occupation numbers n^. 
The occupation numbers are widely used in generalizations of DFT when extending DFT to finite 
temperatures, when considering excitations energies and so on, for a short review see Ref. [29,31]. 
Taking into account n^, we can use the well-known Landau equation of the Fermi liquid theory [32] 

5E , , 

We note that there exist close relations between the Fermi-liquid theory and DFT [11], while the 
Fermi-liquid theory is applicable in principle essentially to all types of Fermi systems including finite 
objects such as atoms, molecules, clusters and fuUerenes. It is difficult to provide an example of a 
Fermi system, which cannot at all be described in the frame of this theory (sec e.g. [32,33]). 

In the functional E\p,n^ both of the parameters, p and nj, are independent variables and deter- 
mine different features of the system under consideration. The ground state energy as a functional of 
the density determines the system's behavior in an external field, while the ground state energy as a 
functional of the occupation numbers characterizes the single-particle spectrum. On the other hand, 
the functional F\p\ given by Eq. (18) depends only on the density. One can use some approximations 
to construct E\p^n^. For instance, writing the density in explicit form given by Eq. (22), we find 
that the functional F\p\ depends upon the occupation numbers [18,19,34] 



F[p, Ui 



(36) 



while £'[p, nj] is given by Eq. (17). In that case, using Eq. (35) one can show that the single-particle 
energies coincide with the eigenvalues £j of Eq. (21) [34]. 

Let us show that the procedure based on Eq. (36) may lead, however, to wrong results. Consider 
a large uniform system of fermions, so that it is possible to introduce the density of states at the 
Fermi-level. In this case, the single-particle potential is a constant, the eigenvalues of Eq. (21) are 
the eigenvalues of noninteracting fermions with bare mass M, and the density of states is equal to 
that of noninteracting fermions. This result is incorrect being related to the fact that the functional 
form i?[p, nj] given by Eq. (36) was artificially chosen. In contrast, as it is seen from Eq. (5), the 
HF ground state energy depends on the occupation numbers quite naturally and leads to the correct 
results when calculating the single-particle spectra [32,35,36]. Inserting Eq. (5) into Eq. (35), we 
can calculate the HF single-particle spectrum ef^ that coincides with the eigenvalues of Eq. (6). 
In case of a large uniform system of electrons, one can check that the density of states tends to 
zero [35,36]. Upon comparing this result with the result of the KS scheme, one may suspect an 
inconsistency in the KS approach. An inconsistency like this was indeed suspected in KS approach 
because Y^iUiei ^ J2i''^iEf^ [18,19], with being calculated in the scheme of [34]. 
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In order to show that the inconsistencies mentioned above are fictitious and to illustrate how 
to calculate the single-particle energies within the DFT, we choose the simplest case when the 
functional is given by Eq. (30). As we shall see, the single-particle energies coincide neither with 
the eigenvalues E^^^ of Eq. (6) nor with £j of Eq. (21). To proceed, we use a method developed in 
[15,28,37,38]. The linear response functions xo, and the density p(r) depend upon the occupation 
numbers. Thus, one can consider the ground state energy as a functional of the density and the 
occupation numbers 



E[p{r),ni] = 7fc[p(r), rii] + FH[p(r), n^] + F^[p(r), n^] + J Kxt(r)p(r)dr. 



(37) 



Here is the functional of the kinetic energy of nonintcracting KS particles. Substituting Eq. (37) 
into Eq. (35) and remembering that the single-particle wave functions and eigenvalues arc given by 
Eq. (21), we see that the single particle spectrum can be represented by the following expression: 



Xo(ri, r2, iw) + 27rp(ri)5(w)(5(ri - ra) 



ri - r2 



dwdvidv2 
2^r ■ 



(38) 



The first and second terms on the right hand side in Eq. (38) are determined by the derivative of 
the functional with respect to the occupation numbers n^. To calculate the derivative we consider 
an auxiliary system of non-interacting particles in a field U{r). The ground state energy Eq of this 
system is given by the following equation 



< = Tfe + y U{v)p{r)dv. 

Varying Eq with respect to the occupation numbers, one gets the desirable result 

8E^ STk 



(39) 



+ < (l)i\U\(j)i >, 



(40) 



provided U = Vh + + v. The third term on the right hand side of Eq. (38) is related to the 
contribution coming from defined by Eq. (27). In the considered simplest case when we take 
into account only T^, Fh and F^ functionals, the infraction enters as a linear factor, with Fh and 
F^ being linear dependent on the strength. If we omit the inter-electron interaction given by Eq. 
(15) we directly get from Eq. (37) = as it must be in the case of nonintcracting system of 
electrons. Note that it is not difficult to include the correlation energy in the simplest local density 
approximation 



Fcip^rii] = j p{r)ec{p{r))dr. 



Here the density p(r) is given by Eq. (22) and the correlation potential is defined as 



Voir) = 



Sp{r) 



(41) 



(42) 



Then Eq. (37) takes the form 
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E[p{r),ni] = Tfe[p(r), n,] + Fj,[p(r), n,] + F,[p(r), n,] + Fe[p(r), n,] + J Kxt(r)p(r)dr. (43) 



Varying £^[p(r),ni] with respect to the occupation numbers rii and after some straightforward calcu- 
lations, we obtain the rather simple expression for the single particle spectrum 

Si = £i- < 0i|T4|0i > - ^rik 

k 

Here we employ Eq. (40) and choose the potential U asU = Vh + Vx + Vc+v to calculate the derivative 
STk/Srii. Approximating the correlation functional Fc[p,ni] by Eq. (41) we simplify the calculations 
a lot, preserving at the same time the asymptotic condition, {V^ + K)r-^oo ~^ —1/r. This condition 
is of crucial importance when calculating the wave functions and eigenvalues of vacant states in the 
frame of the KS approach [28,37]. Note, that these functions and eigenvalues that enter Eq. (44) 
determine the single particle spectrum e^. This spectrum has to be compared with the experimental 
results. 

Single particle levels given by Eq. (44) resemble the ones that arc obtained using HP method. 
If the wave functions (pi were solutions of the HE equations, the correlation potential Vc would be 
omitted, then the energies would exactly coincide with the eigenvalues But this is not the case, 
since (pi and Ei are the eigenfunctions and eigenvalues of Eq. (21). Therefore, the energies coincide 
neither with HE eigenvalues given by Eq. (6), nor with the ones of Eq. (21). As we have seen in 
Section III, in case of the He atom T^^[p\ = Tfc[p] and F^^[p] = F^^[p]. Then it follows from Eq. 
(44) that Y^i^i^i = Hi^iE^^ = Y^i^i^i in accordance with numerical results obtained in [18]. In 
case of other atoms we cannot expect this equality to hold because of the non-locality of the HE 
single-particle potential. This result is in agreement with the numerical calculations [18,19] and has 
nothing to do with contradictions between the HE method and DET. 

It is worth noting that if the exchange functional F.x[p] is treated within the local density ap- 
proximation, as it is done for Fc[p], that is one uses Eq. (36) to construct £'[p, n^], then the second 
and the third terms on the right hand side of Eq. (44) cancel each other, leading to = £j in 
accordance with [34]. Therefore it is of crucial importance to keep the proper representation of the 
functional F,j\p] given by Eq. (27). We also anticipate that Eq. (44) when applied to calculations of 
many-electron systems such as atoms and molecules will produce reasonable results for the energy 
gap separating the occupied and empty states. In case of solids, we expect that the energy gap at 
various high-symmetry points in the Brillouin zone of semiconductors and dielectrics can also be 
reproduced. 

V. CONCLUDING REMARKS 

We have shown that there are no contradictions in DET related to the fact that the HE-like 
calculations of finite electron systems carried out within the Kohn-Sham theory of DET lead to the 
ground state energy, the absolute value of which is smaller than that calculated within the HE method. 
On the other hand, we have demonstrated that the properly constructed HE functional of the density 
produces the ground state energy, which coincides with the energy obtained within HE calculations. 
While inequalities (1) and (2) came from omitting the correlation energy, which is defined by the 
difference between the sum of the kinetic and exchange energies of a system considered within the 



(ri)0i(r2)0^(r2)0fc(ri) 



ri - r2 



dvidv2. 



(44) 
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KS approach and the HF method, respectively. It is the nonlocal exchange potential entering the 
HF equations that generates this correlation energy. We have shown that this functional F^^ of 
the correlation energy, being the specific feature of a finite system, cannot be approximated by the 
local density approximation based on calculations of a homogeneous electron system. Thus, the 
correlation energy of a finite electron system, which has to include F^^[p\. cannot be obtained from 
considerations of uniform electron systems, although this way of constructing Ec[p\ is generally used. 
We can conclude that the HF method has unique peculiarities which can be useful when investigating 
many-electron systems. 

We have clarified the relationship between the eigenvalues of the single-particle KS equations, the 
HF eigenvalues and the real single-particle spectrum. Again, there are no contradictions based on 
the comparison of the single-particle spectra calculated within DFT and the HF method. We have 
presented a simple equation defining the single-particle spectrum of a many-electron system. It is 
worth noting that all the specific features of the HF method mentioned in Section II are dictated 
by the HF nonlocal potential and therefore give no grounds to seek possible contradictions in the 
Kohn-Sham theory which deals with local single-particle potentials and v representable densities. 
As it was mentioned in Section II, the HF method has some drawbacks. Among them there is a 
well-known consequence of the non-locality of the HF single-particle exchange potential, which leads 
to the fact that the HF single-particle potential acting on the unoccupied states falls of exponentially. 
As a result, the HF potential can support very few unoccupied states. Such a poor situation with 
the number of states develops difficulties in treating the excited states. By contrast, the KS theory 
does not suffer from the drawbacks of the HF method. 
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